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Abstract 

. The longwavelength diffusion coefficient of a critical fluid con- 
fined between two parallel plates seperated by a distance L is strongly- 
affected by the finite size. Finite size scaling leads us to expect that 
the vanishing of the diffusion coefficient as for ^ < L, ^ being the 
correlation length, would crossover to L _1 for £ ^> L. We show that 
this is not strictly true. There is a logarthmic scaling violation. We 
construct a Kawasasi like scaling function that connects the thermo- 
dynamic regime to the extreme critical (£ L) regime. 
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One of the earliest and certainly one of the most frequently used scaling 
function in the whole of critical dynamics is the scaling function for the ther- 
mal conductivity / concentration diffusivity in a single component / binary 
fluid near the gas-liquid critical point / critical mixing point. This is the well 
known Kawasaki function which describes the passage of the diffusion coeffi- 
cient (we will use the binary liquid language) from a long-wavelength, finite 
correlation length regime to a finite-wavelength, infinite correlation length 
regime. When the wavelength is very large (k, the wave number — > 0), the 
diffusion coefficient D diverges as £ 4 ~ d = K d ~ 4 , where £ is the correlation 
length; k, the inverse correlation length and d, the dimensionality of space. 
At the critical point (k — 0), this behaviour changes to k d ~ A . In the physical 
(d=3) situation, we have the passage described by the Kawasaki function 
K(x) [1]: 

= (1) 

with 

K(x) = A, [1 + x 2 + (x 3 - -)tan~ l x] (2) 
Ax 2 x 

If k -»• 0, D(k) = and if re -> 0, D(k) = ^-^k. In the above t? is the 
background critical viscosity. The small divergence of the critical viscosity is 
ignored in this paper. A practical approximation to the Kawasaki function 
is to write it as 

D(k,n) = ^L (3) 
with L having the property that 

L — > — as k — > 

K 
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as k 



8 k k k 

(4) 

The first correction being known to be of the form ^ about the k = 
limit, we can use the approximation L(k, k) = {k 2 + k 2 )^ 1 and the diffusion 
coefficient is 

D(fc,K) = -^(fc 2 + «; 2 )3 ( 5 ) 
6nr) 

This is exact for k — > and over the range of practical values of the | ratio, 
is accurate to about 5%. This makes it a very useful approximation. 

In this paper, we want to look at the diffusion coefficient in a finite ge- 
ometry - the fluid contained between two slabs separated by a distance L. 
Recently Koch and Dohm have explored finite size effects on the diffusion co- 
efficient in three dimensional Ising-like systems [17]. Unlike their situation, 
we will show that our system exhibits a small scaling violation. This is a 
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physical model and hence this violation should be experimentally accessible. 
The hydrodynamic limit is taken to hold i.e. we are in the long wavelength 
limit (i.e. k=0). For L » the diffusion coefficient will be governed by 
D(k) = fj^K- In the finite geometry, the leading correction to this result 
was first obtained by Calvo and Ferrell. The result was a bit of a surprise in 
the presence of a logarithmic factor 

H^L) = k- 1 [1- 1 -^], (6) 

This was later obtained using a different technique by one of the present 
authors [16]. The appearance of the logarithmic term prompted us to explore 
the limit £ L. In this limit, we find (as explained below) 

L( K = 0,L) = A ; „_1 

To arrive at the diffusion coefficient, one needs 

ksT L(k, L) 



D{k,L) 



6tt?7o x(«, L) 



where x( K > L) is the susceptibility in the finite geometry. The limiting form 
of x( K ? L) are x{ K i L — > oo) = k~ 2 and x( K ~^ 0, L) = j^, so that 

D(k,K) -> K (l -___), /or « 1 <L 

07T?7o 2 



1 / 1 * -i ^ r 
Tin—, for k > L 

bnr/o L kL 



We combine the two forms to propose the scaling function 



(7) 



MLn j. 1 ^ + §M1 + A)] ; 

6tT7]oC 

k B T 



D(k,K) ~ en Vo ^ + n^ ) l + ^- L ln(.L+£ 



6tt??o£ 



F(kL) 



(8) 



which is the («, L) analogue of the Kawasaki scaling function, given in Eq.(2). 

In the thermodynamic limit (kL — > oo), F(kL) = 1 and rises to \ln-^ 
as — > 0. The change in kL is brought about by varying £ at a fixed 
L. At kL = 1, F(kL) = 1.57, significantly different from its value in the 
thermodynamic limit. This should make the effect observable. To see the 
existence of the logarithmic terms, one would need a fairly high degree of 
accuracy. We now proceed to outline the technical details. 
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The free energy functional that governs the static fluctuation of the order 
parameter ip can be taken to be quadratic and is given by 



F = fd^rf rf4E^ 2 + EW.f 
J0 Z i=l z i=0 

| d D -V[^ 2 (f, z = 0) + ^ 2 (f, z = L)} 



c 
A 



(9) 



Since the anomalous dimension index r\ plays an insignificant role in the 
study of dynamics, we can work with this quadratic expression for F. Here 
the geometry is restricted in the z-direction between z = and z = L. The 
second term on the right hand side of Eq.(9) is a surface contribution, V is 
a constant and A, an extrapolation length. For A — > 0, if) 2 must vanish at 
z = and z = L in order to satisfy the Dirichlet boundary conditions. We 
set A = 0. 

We introduce the Fourier transform of ip(f, z) through the relation 

m z) = JL J z)e^d 2 k (10) 

where k is the wave-vector in the 2-dimensional space. The two-point corre- 
lation function G(k, zi, z 2 ) is the solution of the differential equation 

(^- k2 - k 2 )G{K zi, z 2 ) = -5( Zl - z 2 ) (11) 

with G(k, zi, z 2 ) vanishing at z = and z = L and is given by 

,r» . sinhaz < sinha(L — z>) . , 

11 2 lasinhaL 

with a 2 = k 2 + k 2 . For L — > and ^ > z 2 , Eq.(12) approximates to 

z<(L - g>) 

2L (13) 

At this point it is easy to check that x{ K ^L) = j J G(k = 0, z±, z 2 )dz± dz 2 
gives the results mentioned above Eq.(7). 

Let us now focus our attention on dynamics. For this we need to introduce 
the equation of motion for -0-field. We take this to be a Langevin equation 
where the potential corresponds to the free energy functional of Eq.(9). For 
non-conserved -0-field, this reads as 

^ = -T(k 2 + k 2 - z, t) + N(k, z, t) (14) 

where N(k, z,t) is the noise, characterized by the correlation 

< N(k 1 ,z 1 ,t 1 )N(k 2 , z 2 , t 2 ) >= 2T5(k 1 + k 2 )S{ Zl - z 2 )S(t 1 - t 2 ) (15) 
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If the order parameter field ip is conserved, there will be an additional factor 
(k 2 — multiplying T in Eq.(14) [2]. The dynamic correlation function 

reads as 

C(k,z 1 ,z 2 ) = < ip(k, zi, uj)i/j(—k, z 2 , —oj) > 

= Lj dz'dz"R + ( Zl ,z')R4z 2 ,z") < N(z')N(z") > 
2 r 

= -J dz'R + (z 1 ,z')R-(z 2 ,z') 

(16) 



where 



with 



sinha±z < sinha±(L — z>) 

R± = ; (1' ) 

2a±sinha±L 



a 2 ± = ±f + k 2 + k 2 (18) 

Now, we note that the current j(f, z, t) associated with a transport process 
is in general a bilinear combination of the order parameter field or a com- 
bination of the order parameter field with a secondary field. Then Kubo's 
formula yields the Onsager coefficient corresponding to the current j(r, z,t) 
as 

A oc Yjrf, J < z u ti) ■ j{f 2 , z 2 , t 2 ) > d 2 r l d 2 r 2 dz l dz 2 dt 1 dt 2 (19) 

where V 2 is the volume in the two dimensional space. 

For the binary liquid system, the order parameter field tp is the density 
difference between the liquid and gaseous phase and the relevant current is 
j = ifjv, where v is the velocity field. For the liquid-gas system, this current is 
proportional to the entropy current and so Kubo's formula yields the thermal 
diffusivity. For a binary mixture the current is the mass current and Kubo's 
formula yields the mass diffusivity. Using j = ipv in Eq.(19), we have 

J < ip{r 1 ,z 1 ,t 1 )v(r 1 ,z 1 ,t 1 ) ■ ?p(r 2 ,z 2 ,t 2 )v(r 2 ,z 2 ,t 2 ) 



A oc y-~ I < Zi.tt )^^. z t Ji) ■ ir(ro. ZoJo)r(r-->. zo.to) > 

xd 2 rid 2 r 2 dzidz 2 dtidt 2 



(20) 



Then the decoupled mode approximation enables us to write the correlation 
function in Eq.(20) as a product of two correlation functions [3] viz. 

< ^(fx, z 1 ,ti)tp(f 2 , z 2 ,t 2 ) > and < v(f 1 ,z 1 ,tl) ■v(f 2 ,z 2 ,t 2 ) > and so finally 
we have 

A oc — / C^(k, zi, z 2 ,u>)C vv {—k, z±, z 2 ,uj)d 2 kdujdzidz 2 (21) 

where and C vv respectively stand for the order parameter and velocity 
correlation functions. 
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We now note that the time scales associated with the velocity field and 
the density field are very different. The velocity field relaxes much faster and 
in the time scale (ti — t 2 ), the density field changes hardly. As a result the 
order parameter correlation function C^{f\ — f 2 , z±, z 2 , ti — t 2 ) can be taken 
to remain at its static value and hence one needs the zero frequency limit of 
the velocity correlation function C vv . Therefore Eq.(21) reduces to 

X(X ~lI d2kdz ^ dz 2 C ^ StatiC (ki zu z 2 )C vv (-k, z u z 2 ,uj = 0) (22) 

The static correlation function for the order parameter ip is given by 
Eq.(2). Let us now proceed to evaluate the zero frequency limit of the velocity 

— * 

correlation function c vv (— k, z±, z 2 , uj). The relaxation dynamics of the v- field 
is governed by 

Ov 2 -* -* 

— = ~T v (k 2 - g^)v(k, z, t) + N v (23) 

The velocity field is solenoidal and so V • v = 0. The above constraint tells 
us that we should work with a field A, such that v = V x A. The correlation 
function Caa for the A-field is given by Eq.(12) with 
k = ^ T = T c . The velocity correlation function follows from 

c - = ( * 2 + ^ )c - (24) 

In the limit L — > 0, it is found to be 

a "=2T [ 2W (4Z2 " Zl)] (25) 

Using Eqs.(13) and (26), the Onsager coefficient, in the three spatial 
dimensions, is 

A oc -jr J (fkdzxdz^^ik, z u z 2 )C vv {-k, z u z 2 , uj = 0) 
= -j- f d 2 k f L [f Zl + f L ]dz 2 

L J JO JO Jz-i 

1 , * z 2 (L — Z\) . 



/; J2 

4rL 4 

since 
1 5 rd 2 k fL 



h! d -wlo L{L - Zl)dZl So Z2(AZ2 ~ Zl)dZ2 ' 

/ dzi[ + }dz 2 = 2 dzi dz 2 
Jo Jo Jz{ Jo Jo 



4rL 4 6 

d 2 k 



Jj^Jo 



96rJ k 2 



(26) 
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Therefore, for L — > 

A oc L J - ~ Lin- (27) 

This establishes the point we wanted to make in Eq.(7). Normalizing to the 
same prefactor for the large k and small k limits leads to the coefficients in 
Eq.(7). 
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